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The electron-phonon coupling and the corresponding energy exchange was investigated experi¬ 
mentally and by ab initio theory in non-equilibrium states of the free-electron metal aluminium. 
The temporal evolution of the atomic mean squared displacement in laser-excited thin free-standing 
films was monitored by femtosecond electron diffraction. The electron-phonon coupling strength was 
obtained for a range of electronic and lattice temperatures from density functional theory molecu¬ 
lar dynamics (DFT-MD) simulations. The electron-phonon coupling parameter extracted from the 
experimental data in the framework of a two-temperature model (TTM) deviates significantly from 
the ab initio values. We introduce a non-thermal lattice model (NLM) for describing non-thermal 
phonon distributions as a sum of thermal distributions of the three phonon branches. The contribu¬ 
tions of individual phonon branches to the electron-phonon coupling are considered independently 
and found to be dominated by longitudinal acoustic phonons. Using all material parameters from 
first-principle calculations besides the phonon-phonon coupling strength, the prediction of the en¬ 
ergy transfer from electrons to phonons by the NLM is in excellent agreement with time-resolved 
diffraction data. Our results suggest that the TTM is insufficient for describing the microscopic en¬ 
ergy flow even for simple metals like aluminium and that the determination of the electron-phonon 
coupling constant from time-resolved experiments by means of the TTM leads to incorrect val¬ 
ues. In contrast, the NLM describing transient phonon populations by three parameters appears 
to be a sufficient model for quantitatively describing electron-lattice equilibration in aluminium. 

We discuss the general applicability of the NLM and provide a criterion for the suitability of the 
two-temperature approximation for other metals. 


I. INTRODUCTION 

The interaction between electrons and lattice vibra¬ 
tions is central to both ground state as well as out-of- 
equilibrium properties of solids. The huge phase space of 
elementary energy excitations of electronic and phononic 
origin renders a universally valid description of electron 
and phonon distributions virtually impossible. The most 
common approach is to assume thermal Fermi-Dirac and 
Bose-Einstein distributions for electrons and phonons, 
respectively, and to describe the non-equilibrium quan¬ 
tum statistics by a Boltzmann equation (see e.g. [U, 01). 
These approximations are employed in the description 
of near-equilibrium phenomena like electrical or thermal 
conductivity as well as in the treatment of states with 
pronounced non-equilibrium distributions. 

Non-equilibrium states occur after impulsive perturba¬ 
tion of a material, for instance, after excitation with fem¬ 
tosecond laser pulses. For the case of metals, Anisimov 
et al. introduced an empirical two-temperature model 
(TTM) describing the energy transfer from the excited 
electrons to the lattice with a single electron-phonon 
coupling parameter Si- Allen refined this approach 
by developing a microscopic theory and derived the 


electron-phonon coupling parameter as the spectral in¬ 
tegral of the zero-temperature Eliashberg spectral func¬ 
tion a^F{uj) comprises the phonon den¬ 

sity of states Fiuj) and the phonon frequency-dependent 
electron-phonon coupling which incorporates all 

allowed scattering processes of electrons with phonons 
of frequency w. While Allen’s work provides a link be¬ 
tween the Eliashberg theory of conventional supercon¬ 
ductors and the electron-lattice equilibration in laser- 
excited metals, the underlying assumptions, in partic¬ 
ular the independence of on the electron temperature 
and the assumption of thermal distributions of electrons 
and phonons at all times, are questionable in the case 
of pronounced non-equilibrium situations sa. This is 
rationalized by comparing the phonon density of states 
(PhDOS) and the Eliashberg function, see the results 
of our DFT calculation in Fig. [TJ in a thermal state, 
phonons follow Bose-Einstein (BE) statistics (thin lines, 
plotted for three temperatures), and the phonon occu¬ 
pation is given by the phonon DOS multiplied with the 
BE statistics. The coupling of electrons to phonons is 
given by the Eliashberg function and is strongest for the 
high-energy phonons, whose occupation is the lowest in a 
thermal state. This implies that energy transfer between 
electrons and lattice in non-equilibrium states leads to 
transient non-thermal phonon distributions. 
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FIG. 1. Normalized phonon density of states (PhDOS) and 
Eliashberg function of aluminium, calculated with DFT and 
the abinit package. The solid green lines show Bose-Einstein 
statistics for three lattice temperatures. 


pling of metals quantitatively in equilibrium and non¬ 
equilibrium situations. Among the first were measure¬ 
ments of the phonon lifetime by neutron diffraction Q, 
which, however, depend on electron-phonon and phonon- 
phonon scattering rates. The Eliashberg functions of 
metals have been determined by electron tunneling spec¬ 
troscopy 3 , which have to be done in the superconduct¬ 
ing state at low temperatures. Angle-resolved photoelec¬ 
tron spectroscopy has been used to measure the electron- 
phonon coupling near surfaces through the experimental 
determination of the electron self-energy [13. 

The electron-phonon coupling in non-equilibrium situ¬ 
ations has been studied by time-resolved measurements. 
Subsequent to impulsive excitation of electrons by a fem¬ 
tosecond laser pulse, the electron-phonon coupling gov¬ 
erns the dynamics of equilibration of electrons and lat¬ 
tice. The most common approach is time-domain ther¬ 
moreflectance measurements [HHil. While these mea¬ 
surements are easy to implement, the relation between 
spectroscopic observables, e.g. the reflectivity, and phys¬ 
ical quantities of the non-equilibrium system is, in gen¬ 
eral, nontrivial. In particular, the respective contribu¬ 
tions of electronic excitation, increasing phonon popula¬ 
tion and thermal expansion to the transient reflectivity 
are per se unknown and approximated by linear depen¬ 
dencies. Time-resolved photoelectron spectroscopy pro¬ 
vides a more direct access to transient electron distribu¬ 
tions which can be related to the electron-phonon interac¬ 
tion [l3 |. In the case of simple metals, however, its intrin¬ 
sic surface sensitivity makes a disentanglement of elec¬ 
tron relaxation due to coupling phenomena and ballistic 
electron transport out of the probed sample volume virtu¬ 
ally impossible EH. In contrast, time-resolved diffrac¬ 
tion provides direct information on the vibrational exci¬ 
tation and long-range order of the crystalline structure. 
In fact, the first studies of photo-induced phase transi¬ 


tions with time-resolved electron diffraction have been 
performed on thin aluminium samples [l^ Hsj] . How¬ 
ever, the quantitative investigation of electron-phonon 
coupling in simple metals requires a high temporal reso¬ 
lution not available in time-resolved electron diffraction 
until recently. 

In this article, we report a combined experimental 
and theoretical study of the energy transfer from photo- 
excited electrons to the lattice in the quasi-free-electron 
metal aluminium by femtosecond electron diffraction 
(FED) and density functional theory (DFT). The ex¬ 
perimental and theoretical approaches are described in 
sections El and uni In section IlYl we first apply a two- 
temperature model to extract electron-phonon coupling 
constants from the diffraction data and find a signifi¬ 
cant deviation from the values obtained from ab initio 
calculations. In section El we introduce a non-thermal 
lattice model allowing for an approximate description of 
non-equilibrium phonon distributions. The predictions of 
this model are compared to the experimentally observed 
atomic mean squared displacements (MSD) and a quanti¬ 
tative agreement of DFT calculations and measurements 
is found. 


II. EXPERIMENTAL METHODS 

The femtosecond electron diffraction measurements 
were conducted with a compact electron diffraction setup 
described in detail elsewhere [13 . We use near-infrared 
pump pulses centered at 800 nm with a duration of ap¬ 
proximately 50 fs (FWHM) to photoexcite 30 nm thick 
films of free-standing, polycrystalline aluminium (Plano 
GmbH) and femtosecond electron pulses to probe their 
atomic structure at different delay times. The pump spot 
size is 250 ^m FWHM and about twice as big as the 
probe spot. The energy of the pump pulses was varied 
to achieve absorbed energy densities between 125 J/cm^ 
and 840J/cm^. The repetition rate of the laser was re¬ 
duced from 1 to 0.5 kHz for the highest fluences to allow 
for thermal relaxation back to room temperature between 
laser pulses. The probe electrons have a kinetic energy 
of 100 keV and the pulses contain a few-thousand elec¬ 
trons. According to numerical simulations, the electron 
pulse duration at the sample position is approximately 
100 fs [ 13 . The diffracted electrons are detected with an 
electron camera and the images are angularly integrated 
to retrieve radial averages, shown in Fig.j^Ka) (diffraction 
image in the inset). Comparison to powder x-ray diffrac¬ 
tion data reveals that the polycrystalline films exhibit 
partial texture with preferred (lOO)-direction. 

For each pump-probe delay, we fit an empirical back¬ 
ground function (Lorentzian -|- fourth order polynomial) 
and pseudo-voigt line-profiles to the peaks in the radial 
averages. The fitting of background and Bragg peaks is 
refined iteratively. The temporal evolution of the rela¬ 
tive intensity of three selected Bragg peaks, as well as 
changes in the integrated background between scattering 
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function of 150 fs FWHM to account for the system re¬ 
sponse. This yields a time constant of r = 350 ± 45 fs 
for the dataset shown. The background rises with a 
time constant of r = 270 ± 20 fs. We note that these 
time constants are significantly shorter than previously 
reported values for the lattice heating time constant of 
(550 ± 80 fs, mi), the Bragg peak decay time of alu¬ 
minum (l.OiO.lps, [12) and the rise of the inelastically 
scattered intensity (0.7 ± O.lps, [12) at comparable ex¬ 
citation levels. The solid lines in Fig. [2Ib) are derived 
from the global fit by inversion of Eq. ([5]) and reproduce 
the behaviour of each diffraction peak. 


FIG. 2. (a) Diffraction image of 30 nm thick polycristalline 
aluminium (inset) and radial average, retrieved by angular in¬ 
tegration. Peaks with multiple labels are too close together to 
separate experimentally, (b) Evolution of the relative inten¬ 
sity of three selected Bragg peaks and the background (BG) 
of diffusely scattered electrons with respect to negative time 
delays for an absorbed energy density of 440 J/cm®. The solid 
lines represent a global fit to the data (see text for details). 


vector s = 0.2A“® and 1.95 A“^, are plotted as a function 
of pump-probe delay in Fig. [2Kb) for an absorbed energy 
density of 440 J /cm®. The intensity of all peaks decreases 
after photoexcitation at t = to and reaches a new steady- 
state in less than one picosecond. On the same timescale, 
the number of diffusely scattered electrons increases, re¬ 
sulting in a dynamic background signal. 

The intensity of scattered electrons into a Bragg diffrac¬ 
tion peak depends on the atoms’ movements, i.e. their 
MSP (y?) from their equilibrium position (see e.g. 
Ref. [12) 

7 =/o ■ exp , (1) 

where Iq is the scattered intensity of a perfect lattice 
and s = 2 sin (0)/A denotes the scattering vector. The 
observed decrease in intensity can thus be attributed to 
an energy transfer from excited electrons to the lattice, 
leading to a higher MSP of the atoms. 

The measured relative intensities of the Bragg peaks 
Irei{t) = I{t)/I{t < to) can be represented in a form that 
does not depend on the scattering vector and therefore 
allows for averaging of all fitted peak heights 

= (”“>«) - ("")(* < <.). (2) 

We include the (111), (220) and (311) reflections as well 
as the averaged (331)/(420) and (531)/(442) reflections, 
which are too close together to be separated in our ex¬ 
periment, in the data analysis. As the analysis only 
relies on relative changes of diffraction intensities, the 
partial texture of the samples does not affect the time- 
dependent signal. We fit the quantity of Eq. © with a 
mono-exponential function, convoluted with a Gaussian 


III. THEORY 


The time evolution of a system of electrons and 
phonons including electron-phonon coupling can be stud¬ 
ied theoretically by means of a set of Bloch-Peierls- 
Boltzmann equations as proposed by Allen [2 . This sys¬ 
tem of equations conserves the energy that is transferred 
between electrons and phonons. In the spirit of the two- 
temperature model, the change in temperature is deter¬ 
mined by the energy transfer from electrons to the lattice 
Zep and vice versa. The electron-phonon energy transfer 
rate is then given by a moment of the Bose distribution 
function nB{q,t) for the phonons [2 




frr rr dEph{t) ^ ^ ..dnsiq^t) 

{Te,Ti,t) = —— ^ ^huj{q,Te,Ti,t) - — -. 


( 3 ) 

Note that the phonon energies huj and therefore the en¬ 
ergy transfer rate depend on the temperatures of the 
system. Z^p can be obtained from the appropriate Boltz¬ 
mann equation and reads in terms of the electron-phonon 
matrix element [2 


Zep{T,,Ti,t) = 

qk 

x5{e{k) - e{k') + hhj{q,Te,Tut)). (4) 


The phonon and electron wave vectors are connected via 
k — k' = q and the difference in single-electron energies 
e(fc) has to match the phonon energy huj{q). The thermal 
factor S{k, k') is given by 


S{k, k\ t) = [f{k, t) - f{k', t)] nsiq, t) 

-f{k' ,t)[l- f{k,t)] (5) 


and accounts for Pauli blocking in the scattering process 
of electrons via the Fermi functions f{k, t). The standard 
approach in calculating this expression is to introduce 
the Eliashberg function as a generalized electron-phonon 
matrix element [2 


a'^F{e,e' ,u}) 


2 

HN^giep) 


X “ ‘^(9))^(e - e{k))S{e' 

kk' 


e{k')). (6) 
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Here, giep) denotes the value of the electronic density 
of states (DOS) at the Fermi edge. We observe that 
the Eliashberg function varies on a meV scale compared 
to the usual eV scale for the electronic DOS. Following 
Wang et al. [1^ , we therefore introduce the first approx¬ 
imation 




g(£)g(£0 

gHep) 

g(g)g(gO 

gHer) 


a‘^F{£F,£p,Uj) 

a^Fiu}). 


(7) 


We also use, based on the same argument of scale between 
the phonon energies and the electronic energies 


fe{£ -huj)- /e(e) = -ftw 


dfeje) 

de 


and g{£F + fu^){q)) ^ gisp)- This results in 
2TTNr 


Zep{Te,Ti,t) = - 


■^ep 

CO 


( 8 ) 


fl(ep’) 

X I duj{hio)'^a'^F{uj,Te,Tpt) n%{uj,Te) - n^p{u},Ti) 


OO 

X J deg'^{e) 


dfe{e,Te) 

de 


(9) 


The Bose functions contain the electron temperature in 
rig and the lattice temperature in n^. Input quantities 
are therefore the electron DOS g^e) and the Eliashberg 
function a^Fioj). 

The heat capacities can be obtained from the various 
density of states. For the electrons, one finds 

OO 

Ce(T„t)= J ^Ii^^g(e)ede. (10) 


Similarly, the lattice heat capacity is given by 

OO 

f dnsifu^^Tut) 

= J - — - F{uj)huduj, (11) 

— OO 

where we denote the phonon DOS with F{u}). 

The calculated energy-transfer rates Zep allow to con¬ 
nect the microscopic theory to the diffraction experi¬ 
ments. The standard approach in describing the tem¬ 
poral evolution of the excited material is the two- 
temperature model, which assumes all phonons to fol¬ 
low Bose statistics at all times. In the following sections, 
we compare the TTM to our data and then expand the 
model to investigate the implications of a non-thermal 
lattice on the temporal evolution of electron and phonon 
distributions. We assume the electronic subsystem to be 
thermal at all times in both models, since the electronic 
potential does not significantly change for excitations up 
to 6 eV and the non-equilibrium electron distribution 


induced by the laser pulses relaxes on the timescale of 
10 fs [H, [ 2 ^, which is much shorter than the timescales 
of interest of this work. Additionally, we assume that 
the electronic excitation equilibrates spatially on the sub- 
100 fs time scale due to fast electron transport [l^ in the 
surface normal direction. 


IV. TWO-TEMPERATURE MODEL 

In the TTM, the temporal evolutions of the thermal 
electron and phonon distributions are described by two 
coupled differential equations with the energy flow be¬ 
tween the subsystems being proportional to the electron- 
phonon coupling parameter Gep = Z^pliT^ — Ti) and the 
temperature difference Te{t) — Ti{t): 

BT 

Ce{Te)-^ = Gep ' (Te - Ti) + f {t - to) , (12a) 

Gi{Ti)^=Gep-{Te-Ti). (12b) 

The function f(t — to) models the energy input from the 
laser to the electrons. Its integral is the absorbed en¬ 
ergy density up and its temporal shape is modeled by a 
gaussian function of 50 fs FWHM centered at t = to- 

The predictions of the TTM can be compared to 
the time-resolved diffraction data by extracting the 
atomic MSD from the measured Bragg-peak intensities 
according to Eq. m, and the subsequent conversion 
to lattice temperature using a fourth-order polynomial 
parametrization of the (equilibrium) Debye-Waller B- 
factor [ 2 ^ . A modification of the B-factor due to phonon 
renormalization caused by electronic excitation is ex¬ 
pected to be marginal as the phonon band structure is 
independent on electronic excitation in aluminium [^ . 
The lattice temperatures reached at long delay times are 
in good agreement with the expected temperature rise 
calculated from the laser fluence and the static heat ca¬ 
pacity of aluminium. 

We continue by solving the differential equations m 
numerically to obtain Te{t) and Ti(t) with Ge and Gi 
taken from the DFT calculations as described in sec¬ 
tion ETH With Gep = Gepiup), Up and the zero time 
of the experiment as free parameters, we minimize the 
difference in T/(t) between experiment and TTM. Three 
data sets are exemplified in Fig. [3Ka) together with the 
best fits of the TTM. For all employed excitation den¬ 
sities, we find that the TTM is able to excellently re¬ 
produce the time-dependent data using the fitted cou¬ 
pling parameters shown in Fig. EKb). The increase of 
Gep(up) with Up suggests a significant dependence of 
Gep on Te, as introduced by Lin et al. [^, since larger 
Up results in transiently larger Te- We thus include a Te- 
dependence of the coupling parameter, Gep = Gep{Te), 
in the TTM. We restrict the dependence to being linear, 
Gep{Te) = go + gi ■ Te, in order to keep the complexity 
of the fit to a manageable level. The result of global fit¬ 
ting of all data sets is shown in Fig. [3Kc) as solid blue 
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FIG. 3. (a) Measured lattice temperature (circles) and simulated lattice temperature with a TTM with optimized parameters 
for three different absorbed energy densities ue- In the bottom, the residual of the TTM and the data of the highest excitation 
density is shown, (b) Electron-phonon coupling parameter Gep as a function of absorbed fluence. (c) Electron-phonon coupling 
dependence on electronic temperature obtained by our measurement. The green line shows Gep(T'e) obtained from DFT 
calculations by evaluating Eq. ©. The dashed light green line is the result from Lin et al. @|. 


line with the shaded area indicating the standard devi¬ 
ation. The data is plotted up to Tg = 4000 K, which 
corresponds to the highest Te{t) reached transiently in 
the measurements. While the experimentally determined 
Gep(Te) is in ver y g ood agreement with the theoretical 
results from Ref. [6(], we observe qualitative and quanti¬ 
tative disagreement with Gep(Te) calculated in this work 
(green solid line in Fig. [31(c)): i) the low-excitation limit 
of the experimental Gep is approximately a factor of two 
smaller than the value obtained from ab initio calcula¬ 
tions; ii) the Tg-dependence of the theoretical Gep{Te) 
is much less pronounced compared to the experimental 
results. 

We consider the good agreement between our ejgperi- 
mental results and the theoretical results of Ref. Q to 
be artificial, as we discuss in detail in the appendix. In 
brief, the function Ggp(Tg) reported in does not purely 
result from first-principles calculations but employs an 
input value Ggp(Tg = 300K) obtained from thermore¬ 
flectance measurements employing a TTM analysis [^ . 
The Mreement in the low-Tg onset of Gep compared to 
Ref. Q therefore does not proof the validity of the TTM 
but only indicates consistency of our data with the ex¬ 
perimental results of Ref. |^. Additionally, we asign 
the exact shape of the Tg-dependence of Gep reported 
in Ref. @ to an insufficient number of fc-points result¬ 
ing in differences in the computed electronic structure of 
aluminium. 

Summarizing this section, we observe a significant dis¬ 
agreement between the Ggp(Tg) extracted from time- 
resolved diffraction data in the framework of a TTM and 
the respective function obtained from first principles cal¬ 
culations. In the following section, we introduce an al¬ 
ternative approach that lifts the assumption of incessant 
thermal distribution of the phonons. 


V. NON-THERMAL LATTICE MODEL 

The comparison of Eliashberg function and phonon 
DOS shown in Fig. |T] reveals an enhanced coupling to 
phonons of high frequencies. This implies that laser- 
induced non-equilibrium between electrons and phonons 
subsequently leads to non-thermal phonon distributions, 
i.e., transiently, the phonon distribution is not resem¬ 
bled by the product of a Bose-Einstein function and 
the phonon DOS. It is well established that non-thermal 
phonon populations occur after photo-excitation of ma¬ 
terials with strong coupling of electrons to specific optical 
phonons. Examples include gra phite exhibiting strongly 
coupled optical phonons and crystals susceptible 

to displacive excitation of coherent phonons (DECP) like 
Peierls-distorted crystals [l^, and certain transition 
metal oxides . The energy flow between electrons and 
vibrational degrees have been described by modifications 
of the TTM. For the case of Peierl-dostorted lattices, the 
TTM has been expanded by an explicit description of the 
optical phonon mode linked to DECP [H, . 

In the following, we introduce a minimal model be¬ 
yond the two-temperature approximation allowing for an 
approximate description of non-thermal phonon distri¬ 
butions in a simple metal where no optical phonons are 
present. For three-dimensional crystals with monoatomic 
unit cells, the phonons uniquely group into three acous¬ 
tic branches with distinct dispersion relations purely due 
to symmetry. The partial contributions of the branches 
to the phonon DOS have been studied by static diffrac¬ 
tion techniques and have been found to be well sepa¬ 
rated in frequency [^, [s^. Similarly, the contribution 
of the three phonon branches to the Eliashberg func¬ 
tion can be identified. Dacorogna et al. calculated the 
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FIG. 4. (a) DFT calculation of the Eliashberg function and its contributions from three phonon branches with 

i = T\,T 2 ,L. (b) Lattice heat capacities of the three branches and of the entire lattice in equilibrium as a function of 

temperature, (c) Electron-phonon coupling calculated for the entire lattice and the individual phonon branches as function of 


Te. 


q- and phonon branch-resolved electron-phonon coupling 
in aluminium along high-symmetry directions and con¬ 
jectured that there is no dominant contribution from 
a specific phonon branch [s^. Our DFT calculations 
performed in the full Brillouin zone, however, reveal a 
pronounced phonon branch-dependence of Gep- Figure 
Ha) shows the individual contributions a^Fiiuj) of the 
three phonon branches to the Eliashberg function, cal¬ 
culated as described in the appendix. There is a clear 
spectral separation of the three contributions, in par¬ 
ticular, the high frequency peak dominating a^F can 
completely be attributed to the longitudinal acoustic 
phonon branch {LA). We use Eq. (I5|) with the partial 
Eliashberg functions a^Fi{uj) to calculate the heat ca¬ 
pacities Cp^i and coupling parameters Gep,i of the three 
branches i = TAi,TA 2 ,LA as a function of Te and com¬ 
pare these to the values of the entire phonon spectrum in 
Figs.Hb) & (c). Whereas the partial heat capacities are 
very similar for all T; relevant in this work, the electron- 
phonon coupling differs by up to a factor four between 
the 7142 and LA branch. 

Consequently, we employ the natural grouping of the 
phonons into its three acoustic branches and introduce 
a non-thermal lattice model (NLM) that accounts for 
non-equilibrium between them. Each phonon subset is 
assumed to be thermal and defined by the temperature 
Tp^i. This is justified by the fact that the spectral shape 
of the a^Fi{uj) agrees significantly better with the re¬ 
spective partial phonon density of states Fi{uj) compared 
to the overall functions a^F and F, as shown in Fig. [5] 
in the appendix. Additionally, for a given temperature, 
the Bose-Einstein function varies less over the spectral 
width of any Fj compared to the full phonon spectrum. 
The NLM therefore is a refinement of the TTM. We 
use Eq. dni) to calculate the energy transfer to subsets 
of phonons by replacing the Eliashberg function a^F{u}) 
with a reduced counterpart, a'^Fi{uj), only containing the 
coupling of electrons to the respective phonon branch. 


The energy flow is determined by their temperature dif¬ 
ference and the evolution of all subsystems is given by 
the set of coupled differential equations 

BT ^ 

Ce{Te)-^ =Y,GepATe) ' {Te - Tp^i) 

i—1 

+ f{t-to), (13a) 

dT 

CpATp.i)-^ =GepATe) ■ {Te - Tp^,) 

+ ^pp,ij ■ {Tp,j — Tp^i). (13b) 

where Gep,i is the electron-phonon coupling parameter 
of the i-tti phonon subsystem to the electrons. The total 
energy transfer between electrons and lattice is then sim¬ 
ply given by the sum of all partial energy transfers. We 
introduced a phonon-phonon coupling parameter Gpp^ij 
to account for phonon-phonon scattering processes, lead¬ 
ing to a direct energy-transfer between branches i and j. 
Due to a lack of experimental or theoretical results, we 
choose to use the same phonon-phonon coupling between 
all branches Gpp^ij = Gpp, the magnitude of which we de¬ 
termine by fitting to the experimental data. Note that, 
in the limit of an infinitely large Gpp, the NLM converges 
to the TTM. 

In the situation of a non-thermal lattice, a conversion 
of the measured MSD to a lattice temperature via the 
Debye-Waller B-factor is invalid. Instead, the contribu¬ 
tions of each phonon branch to the MSD can be calcu¬ 
lated from the respective temperature calculated by the 
NLM. The connection of temperature and MSD is given 
by 


OO 



0 


where M is the mass of an atom. By integrating over 
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FIG. 5. (a) Measured mean square displacement (MSD) as a 
function of delay (circles) and calculated evolution predicted 
by the non-thermal lattice model with a phonon-phonon cou¬ 
pling of 3.5 • 10^^Wm“®K“^ (lines). All other material pa¬ 
rameters were taken from DFT calculations. In the bottom, 
the residuals for the highest excitation density are shown, 
(c) Evolution of the temperatures of electrons (blue) and the 
three phonon branches according to the NLM for an excita¬ 
tion density of 680 J/cm®. The inset shows the full range of 
Tbt). 


the partial DOS Fi{uj) and using the respective tem¬ 
perature Tp^i, we obtain the contribution of all three 
phonon branches (ttf) to the total MSD. With the phonon 
branches being orthogonal, this is given by 

(15) 

i=l 

FigureEKa) compares the MSD determined experimen¬ 
tally with the MSD obtained with the optimized solu¬ 
tion of the NLM, where Ce{Te), Cp^i{Tp^i) and, in par¬ 


ticular, Gep,i{Te) are taken from the ab initio calcula¬ 
tion, and Gpp, ue and to are optimized iteratively. The 
only free parameter affecting the dynamics of the MSD is 
Gpp. We find excellent agreement between data and NLM 
using a phonon-phonon coupling of approximately 3.5 • 
10^^WK~^m“^ for all performed measurements. Even 
though literature values for phonon-phonon coupling are 
elusive, it has been concluded from neutron-diffraction 
measurements that electron-phonon and phonon-phonon 
coupling are of similar magnitude , which is consistent 
with our observation. 

Figure[SJb) shows Tg and Tp ^ of the optimized solution 
of the NLM. The different evolution of the temperatures 
of the three branches is pronounced and shows the tran¬ 
sient non-thermal state of the phonons. As expected from 
the previous discussion, Tp^LA rises faster than the tem¬ 
perature of the transverse branches. Tp^EA overshoots 
before the LA phonons eventually equilibrate with the 
two other branches. The timescale of lattice thermaliza- 
tion is sensitive to Gpp and is found to be approximately 
2 ps. 


VI. DISCUSSION AND CONCLUSION 

Whereas the TTM and the NLM can both be applied 
to describe the temporal evolution of the measured MSD 
with similar quality, the two models differ significantly in 
the electron-phonon coupling constants. Gep extracted 
from the data employing the TTM is a factor of two 
smaller than the value obtained from ab initio calcula¬ 
tions, i.e., the integrated electron-phonon coupling of all 
phonon branches. While assuming thermal phonon dis¬ 
tributions, which is a basic approximation of the TTM, is 
inappropriate in non-equilibrium states, the NLM allows 
for an approximate description of non-thermal phonon 
distributions. Our results imply that the determina¬ 
tion of the electron-phonon coupling of metals with time- 
resolved measurements by means of a TTM analysis may 
lead to systematic errors. In contrast, the excellent agree¬ 
ment between data and NLM employing electron-phonon 
coupling constants from theory suggests that the NLM 
represents a minimal model capturing the physics of the 
electron-phonon coupling in aluminium with sufficient ac¬ 
curacy. 

The important practical difference between the models 
is the predicted dynamics of the energy flow between elec¬ 
trons and phonons. Figure IHKa) & (b) visually summa¬ 
rizes the two-temperature and non-thermal lattice model 
and shows the channels of energy flow. The widths of 
the arrows indicate the effective coupling strengths. Fig¬ 
ure He) shows the evolution of the energy content in 
electrons and phonons added by the laser pulse after ex¬ 
citation of the electrons at t = 0. The energy content 
is calculated from the electronic and phononic temper¬ 
atures given by the two models and the respective heat 
capacities. In the NLM, the energies of the three phonon 
branches are summed to obtain the added energy in all 
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FIG. 6. (a) & (b) Schematic representation of the TTM and the NLM. In the latter, the phononic system is subdivided into 
its branches and the couplings are considered separately. The widths of the arrows are proportional to the coupling strengths 
retrieved from theory and experiment, (c) Temporal evolution of the energy content in the electronic (dash-dotted lines) 
and lattice (solid lines) subsystems for an excitation density of 680 J/cm® according to the NLM and the TTM. The energy 
content of the individual phonon branches is shown as dotted lines. The black dashed line indicates 95% of the energy which 
is eventually transferred from electrons to the lattice. 


lattice degrees of freedom. 

It can be seen that the total energy transfer from 
electrons to the lattice proceeds significantly faster in 
the NLM. The dotted line depicts 95% of the energy, 
which is eventually transferred to the lattice. Between 
the two models, the time at which this level is reached 
differs by 40-60% depending on the excitation density. 
While the TTM is perfectly suitable for describing the 
time-resolved diffraction data, as shown in Fig. El^a), it 
returns an ’effective’ electron-phonon coupling constant 
which underestimates the rate of energy flow between 
electrons and phonons approximately by a factor two. 
This has implications wherever the electron-phonon cou¬ 
pling is used to predict the energy flow between the sub¬ 
systems, for instance in non-equilibrium MD simulations, 
the modeling of laser-induced ablation, shock wave and 
plasma generation, as well as in the description of ul¬ 
trafast magnetization dynamics, where the energy flow 
between electrons, lattice and spins is decisive |38l |. In 
addition, the few-lOOfs timescale of energy transfer in 
aluminium implies that photo-induced structural dynam¬ 
ics, e.g. laser-induced phase transitions, may occur on the 
sub-ps timescale purely due to incoherent energy transfer. 
We note that non-thermal electron-distributions might 
need to be considered additionally for materials exhibit¬ 
ing weak electron-electron scattering [^ . 

The refinement of the description of electron-phonon 
coupling including phonon-specificity is not only relevant 
for the evolution of impulsively generated highly non¬ 
equilibrium states [Tsl. [^. Equally important, in quasi¬ 
stationary non-equilibrium situation like the flow of elec¬ 
trical current, energy dissipation from the electrons to 


the lattice is governed by the same coupling phenom¬ 
ena. In other words, the phonon branch-dependence of 
the electron-phonon coupling may also be applied to the 
description of electron transport with a Boltzmann equa¬ 
tion [^ . 

The experimental basis of this work is the time- 
resolved measurement of the atomic mean-square dis¬ 
placement caused by all phonons of aluminium. The 
energy transfer from electrons to phonons is then in¬ 
ferred on the basis of two models, as it is not an observ¬ 
able itself. However, a more detailed view on transient 
phonon distributions could be obtained by time- and 
momentum-resolved measurements of inelastically scat¬ 
tered electrons [4l| or x-rays from single-crystalline 
thin films of aluminium. Such studies can provide direct 
experimental evidence of transient non-thermal phonon 
distributions, as very recently demonstrated for the nobel 
metal gold (4ll| . 

Heterogeneous coupling of electrons to different types 
of phonons is well established for materials exhibiting 
strongly coupled optical phonons . Our work indi¬ 

cates that phonon-specific coupling has to be considered 
even in simple crystals lacking optical phonons like the 
free-electron metal aluminium. 

The concept of subdividing the lattice degrees of free¬ 
dom into individual phonon branches for describing the 
energy flow in a material, as introduced here for the case 
of aluminium, may be applied to a range of materials. 
In the case of metals, a pronounced spectral dependence 
of suggests that this approach is required, and we 

expect partitioning into phonon branches to be a gen¬ 
erally suitable approach. More complex materials may 
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require partitioning into a larger number of phononic 
subsystems. Ab initio calculations can guide in identi¬ 
fying the appropriate subsystems and provide the partial 
coupling constants, as demonstrated here. 
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Appendix: Computational details 

All simulations were performed with the ab initio sim¬ 
ulation package abinit 13 . The electron-phonon 
cou plin g calculations capabilities are described in Refs. 

B l46l | and are implemented similarly to the one in Ref. 

. We used a Troullier-Martins type norm conserving 
pseudopotential from the fhi database with the outer¬ 
most three electrons taken into account explicitly and 
the remaining 10 electrons frozen in the core [3l- The 
planewave cutoff was Ecut = 13 Ha and the maximum 
angular quantum number is 1 = 2. We checked the re¬ 
sults for the electronic DOS and the Eliashberg function 
with a calculation featuring a harder pseudopotential re¬ 
quiring a plane wave cutoff of 40 Ha and confirmed our 
results. The exchange correlation potential was taken 
in the generalized gradient approximation (PBE-GGA) 
(3 - Additionally, we performed a calculation in local 
density approximation (LDA), which gave a slightly dif¬ 
ferent shape of the Eliashberg function in phonon fre¬ 
quency space but left the first moment of the Eliashberg 
function unchanged. We took the experimental lattice 
constant for fee aluminium a = 7.652 as. The k-point 
grid was unshifted featuring 32 x 32 x 32 points. The 
q-point grid was commensurable with the k-point grid at 
8 x8x8 points. The occupation of the eigenvalues was 
determined by a Fermi distribution with 0.03 eV smear¬ 
ing. 

Such a calculation was performed for the unit cell 
of fee aluminium in order to obtain the values for the 
Eliashberg function. Similarly, the electronic DOS was 
obtained using the tetrahedron method featuring a k- 
point grid up to 128 x 128 x 128 for the ideal lattice at 
Ti = OK. Electronic temperature effects on the DOS and 
the Eliashberg function up to Tg = 20000 K have been 
checked and found to be small. In order to analyze lat¬ 
tice temperature effects on the electron-phonon energy 
transfer, we prepared density functional molecular dy¬ 
namics (DFT-MD) runs with Ni = 32 aluminium atoms 
in a super cell at normal density of p = 2.7g/cc. The 
MD part used a Verlet algorithm and Nose thermostat 



E-Ep [eV] 

FIG. 7. DFT calculation of the electronic density of states in 
aluminium. The dashed line illustrates the ^/E behaviour of 
a free e-gas. 



Energy [meV] 


FIG. 8. Gomparison of the spectral shapes of the partial 
Eliashberg functions a^Fi and partial phonon density of states 
Fi. Each function is normalized to its integral. 


with a time step of At = 0.2 fs. We allowed the ions in 
the supercell to equilibrate over several thousand time 
steps at Ti = 300 K, Ti = 500 K, Ti = 700 K, Ti = 900 K, 
and Ti = HOOK. From the equilibrated part of the runs 
representing fee lattices with irregular harmonic and an- 
harmonic lattice distortions, we randomly took 10 snap¬ 
shots each and subjected these snapshots to the same 
electron-phonon calculations as the initial unit cell. As 
usual, three perturbation directions per atom have to be 
considered in order to build the electron-phonon matrix 
elements, requiring 96 separate calculations per tempera¬ 
ture and snapshot. In addition, due to the size of super¬ 
cell and constraints in the available computer memory, 
we chose to reduce the k-point grid for these cases to 
4x4x4 and the q-point grid to 2 x 2 x 2. We also ran 
simulations with 8x8x8 k-points to check on the error 
introduced this way and found only very minor devia¬ 
tions in the Eliashberg function. As the electronic DOS 
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ill aluminium is already close to the one for an ideal free 
electron gas, we observed only small changes in the DOS 
with an increase in lattice temperature. Using the elec¬ 
tronic DOS we can also obtain the electronic heat capac¬ 
ity needed in the TTM The lattice heat capacity is 
obtained using the same phonons as used for the calcu¬ 
lation of the electron-phonon coupling. 

The partial Eliasberg functions for the NLM can be 
obtained from the phonon line-widths and the partial 
phonon DOS as computed with abinit. A practical 
and simple approximation is given by dividing the total 
Eliashberg function by the total Phonon DOS and using 
the quotient together with the partial phonon DOS to 
obtain partial Eliashberg functions. In Fig. |8] the partial 
phonon DOS Fi and the so-obtained partial Eliashberg 
functions a'^F{uj), normalized by their respective inte¬ 
gral, are plotted. 

The difference in Gep between our result and the one 
from Lin et al. see Fig. [31 is caused by two different 
reasons. Firstly, there is the difference in the electronic 
DOS. The pronounced oscillations in the DOS in Ref. 
in particular near the Fermi level, cause the strong rise 


of Gep(Te) for Te increasing from 0 to 5000 K. As such 
features are absent in the DOS of this work, which is 
in v ery good agreement with other ab initio calculations 
[^[^, the Tg-dependence of is weak. We note that 
the DOS of aluminium of Ref. @ is reproduced by the 
tetrahedron method with 11x11x11 k-points. Sec¬ 
ondly, and more importantly, there is a difference in the 
value for the moment of the Eliashberg function A(a;^). 
Our result A{a;2) = 4.81 ■ 10"^ Ha^, A = 0.47 is ob¬ 
tained entirely on the basis of our DFT calculations and 
agrees well with other DFT calculations (A = 0.43 M), 
(A(a;^) = 4.79 ■ 10 ^ Ha^, A = 0.44) and experimental re¬ 
sults (A^^) = 5.1-10“^ Ha^, A = 0.42) given by Savrasov 
et al. (47|. In addition, our prediction for the supercon¬ 
ductivity temperature as obtained from the values for 
the electron-phonon coupling via the McMillan formula 
[ 5 ^ is Tc = 1.35 K compared to the experimental value 
of Tc = 1.2 K. The value of A(a;^) = 2.51 • 10“^ Ha^ used 
by Lin et al. is not obtained from DFT but is based on 
experimental results employing a TTM in the data anal¬ 
ysis [ 2 ^ and uses the approximation (w^) = with 

9d being the experimental Debye temperature. 
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